Magnetization of mesoscopic disordered networks 
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We study the magnetic response of mesoscopic metallic isolated networks. We calculate the average 
and typical magnetizations in the diffusive regime for non-interacting electrons or in the first order 
Hartree-Fock approximation. These quantities are related to the return probability for a diffusive 
particle on the corresponding network. By solution of the diffusion equation on various types of 
networks, including a ring with arms or an infinite square network, we deduce the corresponding 
magnetizations. In the case of an infinite network, the Hartree-Fock average magnetization stays 
finite in the thermodynamic limit. 
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The problem of persistent currents in mesoscopic rings 
has been stimulated by a few key experiments in 
the recent years. Two types of measurements have been 
deviced, single ring experiments [||j|J^ and many rings 
experiments [0,||J|]. In the second case, the measured 
quantity is an average magnetization (M) while the first 
type of experiment can only give the magnetization cor- 
responding to a given disorder configuration. In the last 
case, the width Mtyp of the magnetization distribution is 
also of interest: M^^^ = (M^) - (M)'^. 

Essentially two types of methods have been used: i) 
analytical methods where the diffusive electronic mo- 
tion is treated in a perturbative way, leading to the fa- 
mous Cooperon diagrams; non-interacting electron the- 
, 12, p^ , p^ , p5[ or Hartree-Fock approximation 
,19 1 have been considered; ii) strictly ID 
models or numerical methods in which either there is 
no diffusive motion or the system size is too small to 
give quantitative results |Q. Up to now, the only re- 
sults for the diffusive regime are given by the pertur- 
bative method. The status of the comparison with ex- 
periments is not completely clear yet but its seems to 
be a reasonable agreement with preliminary recent data 
the typical current being properly well described by 
the non-interacting theory [|9|, |l5| , pT| and the average cur- 
rent being described by the Hartree-Fock approximation 

m. 




a|7 18 



We propose that a new way to get insight into the 
problem is to study other geometries than simple rings. 
In this letter, we calculate analytically the typical and 
average magnetizations of various types of networks, fol- 
lowing the method (i). To do so, we use a semi-classical 
picture to relate the quantities of interest to the return 
probability of a classically diffusive particle. Then, this 
return probability is calculated on different type of net- 
works, giving access to the magnetization. As examples, 



we treat the case of an isolated ring connected to one or 
two arms, and the case of an infinite square lattice. Sev- 
eral experiments are proposed (otherwise specified, h = 1 
throughout the paper). 

In the absence of e— e interactions, a finite contribution 
to the average magnetization comes from the fact that the 
number N of particles is fixed in each subsystem of the 
ensemble . It turns out that this contribution is by far 
smaller than the experimental results. However we will 
discuss it mainly for pedagogial purpose and comparison 
with other contributions. With this constraint on N, the 
" canonical" magnetization is given by [EJ : 



{M^(H)) = -^^{^N^i^,)), 



(1) 



where A is the mean level spacing and {SN'^{fi)) is the 
sample to sample fluctuation of the number of single- 
particle states below the Fermi energy /i. It is an integral 
of the two-point correlation function of the density of 
states (DOS) K{ei — £2) = {p{ei)p{e2)) — Po ■ Po is the 
average DOS. K{e) has been calculated by Altshuler and 
Shklovskii and later in the presence of a magnetic flux 
]lO| , pd 12 1 . A very useful semiclassical picture has been 
presented by Argaman et at, which relates the Fourier 
transform K{t) of K{e) to the classical return probability 
p(r,r,t) for a diffusive particle |l3|. 
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where P{t) = J p{r, r, t)dr. This return probability has 
two components, the purely classical one and the interfer- 
ence term which results from interferences between pairs 
of time-reversed trajectories. In the diagrammatic pic- 
ture, they are related to the diffuson and Cooperon di- 
agrams. The interference term is field dependent and is 
solution of the diffusion equation: 
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[| - ^(V + r', t) = - r') (3) 

From eqs. (0,|^), the average canonical magnetization 
can be related to the field dependent part of the return 
probability: 



{Mn{H)) 
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dt. 



(4) 



Note that the field dependent part of this integral con- 
verges at small times. At large times, the return proba- 
bility is exponentially cut-off as e""*"* where 7 = hD/L^ 
is the inelastic scattering rate. 

Due to the e — e interactions, a larger contribution to 
the average magnetization exists, which has been calcu- 
lated by Ambegaokar and Eckern |16|, in the Hartree- 
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Fock approximation. It can be written as |16 
(Af,e(ff))=-j^ j{n{vf)dv 



(5) 



Where U is an effective screened interaction and n(r) is 
the local density. The integrand is related to the fluctu- 
ations of the local DOS which in turn can be related to 
the return probability One gets: 
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In a similar way, the typical magnetization can also be 
straightforwardly written in terms of K{£) ||l5| , pT|| . By 
Fourier transform, one has: 
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P"{t,H 
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where P"(i, H)\^ ^ d^P/dH\ - d^P/dH^ln- 

To be complete, we remind that the weak-localization 
correction to the conductance of a connected mesocopic 
sample can be also be related to the return probability 
p5|,p6|: 



Aa{r) = {-2/npo)aoC^ir,r) 



(8) 



ctq is the Drude conductivity. The Cooperon C^{r,r, H) 
is the time integrated field-dependent return probability: 



C^{r,r,H)= / p{r,r,t,H)dt 
Jo 



(9) 



It appears that all the quantities of interest are obtained 
as time integrals of the return probability with various 
power-law weighting functions. Noting that P{t) has the 
form Po(i)e^^* and that 



J ^^e-^'dt ^ J°° dj J C^(r,r,H)dr 



(10) 



the different magnetizations can be given in terms of the 
successive integrals of C^(r, r, H): 



(^W))=-^^Jc^^^\r,r,H)dr 
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Kp(H) = ^^, f (r,r,i7)|> 



(11) 
(12) 
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where C- 
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d-fn- d-fi d-y'Cy . These are 
the key equations of this paper since the different magne- 
tizations can be calculated from the knowledge of the re- 
turn probability C^{r, r, H) on the different lattices con- 
sidered and can be deduced from each other or related to 
weak-localization correction by H— or 7— derivatives or 
integrations. 

For the case of weak-localization correction, an exten- 
sive study of this quantity on various lattices has been 
carried out by Dougot and Rammal [ p6| . Considering 
networks made of quasi- ID wires so that the diffusion can 
be considered as one-dimensional, the Cooperon C^(r, r') 
obeys the diffusion equation 



with the continuity equations 

^ . d 2eA „ . ,s « . 



(14) 



(15) 



r, r' are linear coordinates on the network and a, (3 are 
nodes. The term 7 in eq. |lj describes the inelastic scat- 
tering. Integration of the differential equation (|lj) with 
the boundary conditions (|l^) leads to the so-called net- 
work equations which relate C^(a, r') to the neighboring 
nodes l3 M. 
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laf3 is the length of the link (a/3) and jaf3 = 
(47r/^o) Ja -^dl is the circulation of the vector potential 
along this link. Then C~^{r',r') is calculated in terms of 
the C-y{a,r'). Finally, spatial integration give access to 
the magnetizations. 

As an example, we have first considered the case of a 
ring of perimeter L connected to an arm of length b. Such 
a geometry has been considered in the stricly ID case 
without disorder |p7| , p8| , p9| . It is expected that since the 
electrons will spend some time in the arm where there 
are not sensitive to the flux, the persistent current will 
be decreased. From eqs. p^ , |l5| , the function C\{r, r) can 
be straightforwardly calculated on the arm and on the 
ring. Spatial integration gives: 



(Mn) = 7 j- 

TT(po tanh sinh 



sin A'KLp 
+ cosh 4^ 



- cos A'KLp 
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where if = 4>/(t>o. (f> is the flux through the ring, (j)o = h/e 
is the flux quantum and S is the area of the ring. Writing 
this magnetization as (Afjv) = ^^(A/jv)m sin47rm(/J, the 
harmonics content is given by 



2AS — ni arg cosh [cosh ■ 
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Well-known results are recovered when 6 = 0. In the 



case where L ; 
written as: 



L^, the harmonics content can be simply 
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so that in the limit 6 — > oo, the harmonics are reduced 
by a ratio (2/3)'", compared to the ring. However the 
arm has another much more dramatic effect which is to 
decrease the interlevel spacing: A(6) — A(0)L/(L + b). 
This reduction does not exist for (Mee) and Mtyp which 
can also be simply calculated by successive integrations 
on 7. Here we give the result for an arm of length b 3> L^: 



:(00)>m = (3)' 



■(Mee(O)), 



Mtyp{oo)r 



(18) 



(19) 



where (Mee(0))„i are the magnetizations of the isolated 
loop. Interestingly, it is seen that the magnetizations 
(Mee) and Mtyp do not decrease to when 6 ^ 00 but 
saturate to finite values with respective reductions of the 
first harmonics in the ratios 2/3 et 1/2/3. In the limit 
b ^ L^p^ the magnetization should be unchanged if a 
reservoir is attached to the arm ||2^. The case of a ring 
connected to two arms of length b shown in figure 1 can be 
treated in a very similar way. In this case we find that the 
lowest harmonics of the e — e average current is reduced 
in a ratio 4/9 and the typical current is reduced by a 
factor 2/3. We propose that single ring experiments with 
appropriately designed arms could be able to measure 
these reductions. 





FIG. 1. Two geometries considered in this letter 



We now turn to the case of an infinite square lattice 
whose magnetization will be compared with the one of an 
array of isolated rings. The eigenvalues of the diffusion 
equation can be calculated for a rational flux per plaque- 
tte Lp = 4>/'t'a = p/29. Denoting by a the lattice param- 



(16) eter, 77 = a/L^ and (j) — Ha the fiux per plaquette, we 



find that the canonical magnetization per plaquette is 



A d 
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^{ln(4cosh77-e,(0,Ai))} (20) 



where {...} If 
tions of the determinental equation 



ei{9,fi) are the solu- 



detM = det 



Ml 1 ... e*'^ 
1 M2 1 ... 







1 M„ 



(21) 



where M„ = 2cos(4n7riy9 + d/q) — e. M is a matrix asso- 
ciated to the Harper equation known to be also relevant 
for other related problems like tight-binding electrons in 
a magnetic field or superconducting networks in a 

field m. 

The magnetization per plaquette can be compared to 
the magnetization of a square ring of perimeter L = 4a: 

A d 

(Ma,) = ^ ^ ln(cosh 477 - cos Ainp) (22) 

Since A — > for the infinite network, this canonical mag- 
netization density vanishes for an infinite network as it 
was already noticed for a chain of connected rings 
On the other hand, the e — e contribution stays finite in 
the thermodynamic limit. It is given by: 



(Me, 



= Upa^rir/. I {ln(4cosh77-£,;(6',^))}r7d77 

(23) 



and can be compared with the ring magnetization which 
can be cast in the form: 
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(This integral can be calculated explicitely in terms of the 
Lobatchevsky function and it has the Fourier decompo- 
sition found by Ambegaokar and Eckern p6[). Contrary 
to the canonical magnetization, the e — e magnetization 
is an extensive quantity. This magnetization density is 
plotted on figure 2 for the ring and the infinite lattice. 
It is first seen that the network magnetization is contin- 
uous. Although the field dependence of the eigenvalues 
of the Harper equation has a very complicated discon- 
tinuous behavior ( the so-called Hofstadter spectrum). 
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the sum on the eigenvalues has a smooth behavior [ p3[ . 
(Mee) can be calculated easily for large q. In this case the 
dispersion ei{9, fi) is very small and the density of states 
can be approximated by a sum of S functions [p4| . A very 



good approximation of the sum (23) can be obtained by 



replacing the integrals on 9 and ^ by the value of each 
term in the sum taken at Si^ir, tt). 
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FIG. 2. Average magnetization (Mee) of a single ring 
(full lines) and magnetization density of the infinite network 
(dashed lines), for — oo, 4a and a 

It is seen on figure 2 that the network magnetization 
density is about 25 five times smaller than the ring mag- 
netization. Considering that on the array of square rings 
already considered experimentally |Q, the distance be- 
tween rings is equal to the size of the squares, the num- 
ber of squares is four times larger when they are con- 
nected. One then expect only a factor of order 6 between 
the magnetization of the array of disconnected rings and 
the lattice. The width of the magnetization distribution 
scales as \/^fS, S being the area of the network. 

Our results have been obtained in the Hartree-Fock ap- 
proximation and should be corrected by higher order con- 
tributions 1 17 3^. But we do not see any reason why the 
ratio between these two magnetizations could be drasti- 
cally changed. Thus we propose that the measurement 
of the magnetization of a mesoscopic network should give 
similar results to the one of an array of disconnected 
rings. 
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